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Quantum fluctuations of the color flux tube 
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The quantum fluctuations of the flux tube joining two static sources in the confining phase of a lattice gauge 
theory are described by an effective string theory. The predictions of the latter for ratios of Wilson loops of equal 
perimeter do not contain any free parameters, and can be computed exactly for large Wilson loops. We compare 
these predictions with numerical results in 3D Z2 gauge theory, finding complete agreement. 



1. INTRODUCTION 

The flux-tube picture of confinement is now a 
20 year old conjecture, which we can summarize 
in two statements: 

• Two static sources in the confined phase of 
a lattice gauge theory are joined by a color 
fl-ux tube which is responsible for confine- 
ment. 



• The flux tube can fluctuate around its equi- 
librium position; in the rough phase, these 
fluctuations are massless 0. 

These two hypotheses not only provide an intu- 
itive, physical picture of conflnement, but have 
precise quantitative consequences. In principle, 
these predictions can be compared with the ex- 
perimental data for the spectrum of bound states 
of heavy quarks; in practice, precise tests are cur- 
rently possible using numerical simulations. 
The aim of this work is to show that the preci- 
sion now available in Monte Carlo data for the 
simplest lattice gauge theory, namely the 3D Z2 
gauge model, allows one to say a final word about 
the correctness of the fiuctuating flux tube pic- 
ture. 



2. THE PREDICTIONS OF THE FLUC- 
TUATING FLUX TUBE MODEL 

The natural quantity to be studied in this con- 
text is the Wilson loop W{R,T). Suppose two 
static sources are pulled apart at a distance R, 
kept apart for a (euclidean) time T, then put back 
together. If we assume that the flux tube keeps 
its equilibrium position, a straight line joining the 
two sources, the expectation value of the Wilson 
loop will follow the area law: 

{W{R,T)) oce-"^'^ (1) 

where a is the string tension. For the interquark 
potential this corresponds to a linear conflning 
potential: 

V{R) = aR (2) 

If, instead, we allow the flux tube to fluctuate 
around its equilibrium position, we obtain cor- 
rections to the area law that can be compared 
to the results of Monte Carlo simulations. Many 
different models can be written down to describe 
the dynamics of the flux tube fluctuations; how- 
ever, a vast class of these reduce to an effective 
free string model in the infrared region, i.e. for 
large Wilson loop. In this region the expectation 
value of the Wilson loop is given by 
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{W{R,T)) cx e 
Imexpj-I J <f^X\-d^)X^ 



(3) 



*Work done in collaboration with M. Caselle, R. Fiore, F. 
Gliozzi and M. Hasenbusch 



The fields X', (i = 2, . . . , d — 1) are defined on 
the rectangle (0, R) x (0, T) and satisfy Dirichlet 
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boundary conditions. They describe the fluctu- 
ations of the flux tube in the d — 2 directions 
transverse to the Wilson loop. 
The functional integral is gaussian and can be 
performed exactly, using for example C-function 
regularization. The result is ||] 



{W{R,T)) cx e' 
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(4) 



where t — and 77 (r) is the Dedekind function, 
which can be expressed as an infinite product 



(5) 
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in terms of the variable q = e'^^'^'^ . Eq. is 
the prediction of the effective string model for the 
expectation value of large Wilson loops. It should 
be noted that that Eq. (Q) does not contain any 
new adjustable parameters with respect to the 
simple area law (m . 



Figure 1. Comparison of the predictions of the 
free string model with Monte Carlo data for the 
ratios R{L,n). The solid line is the prediction of 
the string model, Eq.(^. 



3. RATIOS OF WILSON LOOPS 



Our goal is to compare Eq.(Q) with the results 
of Monte Carlo simulations of the 3D Z2 gauge 
model. The choice of the model is dictated by 
the need to have very high precision data to ob- 
tain an unambiguous result about the validity of 
the string model. This test could be performed 
by fitting a set of Monte Carlo data for the Wil- 
son loop with Eq.(^). Taking into account the 
presence of a perimeter term, this would be a fit 
with three adjustable parameters. 
However, to obtain a completely unambiguous 
test of our ansatz, we decided to follow a different 
strategy, namely to define an observable quantity 
for which the predictions of the string model do 
not contain any adjustable parameters. We con- 
sidered the following ratios of Wilson loops: 



R{L,n) 



{W{L + n,L~n)) _,„2 
{W{L,L)) ' 



(6) 



Using Eq.(Q) we obtain the prediction of the 
effective free string model for the ratios R{L,n), 
which turns out to depend only on the asymmetry 



ratio t = n/L: 
i?(L,n) = F{t) = 



V 



(■«) 



1/2 



(7) 



4. COMPARISON TO THE MONTE 
CARLO DATA 

We performed a Monte Carlo simulation of the 
3D Z2 gauge model at several values of the cou- 
pling /3, all located in the rough phase and close 
enough to the deconfinement point to be well in- 
side the scaling region. From these simulations 
we extracted the expectation values of the Wil- 
son loops. To evaluate R{L,n), one needs also 
to know the value of the string tension a. How- 
ever, precise evaluations of the interface tension 
in the 3D spin Ising model in the scaling region 
are available. Using the duality between the Z2 
gauge model and the Ising spin model in 3D, we 
can simply plug the Ising model values in the ex- 
pression (1). In this way we use values of a that 
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Table 1 



Monte Carlo 


results for the ratios n) 








n/L 


L 




a 


R{L,n) 


F{n/L) 


0.2 


15 


0.75202 


0.01023 


1.0104(17) 


1.00881 


0.25 


20 


0.75632 


0.004779 


1.0166(10) 


1.01453 


0.33 


12 


0.74883 


0.014728 


1.02881(30) 


1.02901 


0.375 
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0.75245 


0.009418 


1.0403(31) 


1.03940 


0.45 


20 


0.75632 


0.004779 


1.0684(23) 


1.06588 


0.5 


20 


0.75632 


0.004779 


1.0911(27) 


1.09153 
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0.004779 


1.165(11) 


1.17667 
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Figure 2. Finite-size effects for small Wilson 
loop. The prediction of the free string model is 
R{L,L/2) = 1.09153 . . . (straight line). 



are independent from our Monte Carlo samples. 

The results of the comparison are shown in 
Fig.(P and Tab.(^. For each value of the asym- 
metry ratio t = n/L we have displayed the value 
of the ratio R{L, n) with the largest available 
physical size aL^. The solid line is the predic- 
tion, Eq. (^), of the effective free string model. 

Two conclusions can be drawn from these data: 

• The fluctuations of the flux tube are quan- 
titatively relevant: it is easy to see that ne- 



glecting them the ratios R{L,n) would be 
predicted to be all equal to one 

• The free string model, Eq. (||) , describes the 
fluctuations with great accuracy for large 
Wilson loops 

For Wilson loops smaller than a threshold size 
of order aL^ ~ 1, rather significant finite-size ef- 
fects appear. In Fig. (^ we display various data 
for ratios R{L, L/2), i.e. with the same asymme- 
try ratio t = 1/2, and different sizes aL^. 

In conclusion, we have shown that, at large dis- 
tances, the physics of the flux tube fluctuations 
is accurately described by a free string model. 
For shorter distances, finite-size effects are rather 
significant. Precise data for large Wilson loops 
are needed to distinguish between the free string 
model described here and other possible models. 
For example, in Ref.|^ a fermionic string model 
was shown to describe the flux-tube fluctuations 
with an accuracy comparable to or even greater 
than the free, bosonic string model. The high pre- 
cision data now available for large Wilson loops 
allow us to select the bosonic model as the most 
accurate. 
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